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Abstract 

Boson mapping techniques are developed to describe valence bond correla- 
tions in quantum spin chains. Applying the method to the alternating bond 
hamiltonian for a generic spin chain, we derive an analytic expression for the 
transition points which gives perfect agreement with existing Density Matrix 
Renormalization Group (DMRG) and Quantum Monte Carlo (QMC) calcu- 
lations. 

PACS numbers: 75.10.Jm 



Antiferromagnetic spin chains have been the subject of intense interest in recent years, 
largely due to the conjecture by Haldane |I] that chains built from integer spins should 
exhibit a gap in their energy spectrum. The existence of this gap, originally predicted on 
the basis of simple field-theoretic considerations, was subsequently confirmed experimentally 

§• 

Much insight into the properties of antiferromagnetic spin chains has been provided 
by simple models. The field-theoretic non-linear sigma model (NLSM) for example, 

provided the original motivation for Haldane's conjecture. Equally illuminating insight into 
the diverse properties of spin chain systems has been provided by the introduction of the 
Valence Bonds Solid (VBS) state f|. 

Detailed descriptions of these systems, however, have depended on very complex numer- 
ical analyses, using either Quantum Monte Carlo (QMC) simulations || or Density-Matrix 
Renormalization Group (DMRG) methods ||. These "exact" treatments provide striking 
confirmation of the features conjectured by Haldane. 

Ideally, it would be nice to have a simpler method for a reliable quantitative treatment of 
quantum spin chains. Valence bonds provide a physically-motivated starting point for such a 
description. The VBS is the exact ground state of specific spin-chain hamiltonians involving 
quadratic and quartic terms, suggesting that wave functions constructed in terms of valence 
bonds might be good trial states more generally. Unfortunately, such wave functions are 
still too complex, for reasons to be discussed later, to be useful in a variational analysis. In 
the present work, we propose a method that permits valence bonds to be used efficiently in 
variational calculations. We also describe the application of this method to study the phase 
transitions in spin chains governed by the alternating bond hamiltonian 0. 

A useful starting point for the introduction of valence bonds is through the Schwinger 
boson realization of the spin algebra || . The basic idea is to introduce a set of boson creation 
and annihilation operators jj ff and 7i j(T , respectively. These operators create and annihilate, 
respectively, a spin-^ boson with spin projection a = +~ (denoted +) or a = (denoted 
-) at site i. A spin-S system would then involve 25 Schwinger bosons on each site. 
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Typical spin-chain hamiltonians involve spin-spin interactions between nearest neigh- 
bors. It is natural therefore to consider states built up in terms of bonds reflecting these 
correlations. This is the basic idea behind the introduction of valence bonds. When deal- 
ing with an antiferromagnetic spin chain, the key correlations involve nearest neighbors in 
spin-singlet states, which can be represented by the singlet bond 

r I = -^{it+ll+i - - 7i,-7jfi,+) • (1) 

In terms of these singlet bonds, the VBS ground state for a spin-S chain (S an integer) is 
given by 



VBS > = II ( r l) 5 l > • ( 2 ) 

i=l,N 



The state \VBS > is the exact ground eigenstate of the hamiltonian 

H = ^2 Si" ^i+l + q ' ^+l) 2 (3) 

i L 6 J 

involving quadratic and quartic spin operators. For a general hamiltonian, however, it is 
not an eigenstate. Furthermore, when not an exact eigenstate, it is not an especially useful 
trial state. The reason is that the singlet bonds r] from which it is built are not bosons, 
i.e., the corresponding operators do not satisfy boson commutation relations. 

This is a familiar problem in many branches of physics. A standard approach to problems 
of this type, involving dominant pair correlations, is to implement a boson mapping |J. 

The general idea of a boson mapping is to replace the original problem involving pair 
degrees of freedom (in this case, singlet bonds) and the true hamiltonian of the system by an 
equivalent problem involving real bosons and an appropriate effective hamiltonian for these 
bosons. All of the exchange effects between the constituents is transferred to the effective 
hamiltonian, in a mathematically rigorous way guaranteed to preserve the physics of the 
original problem. 

In more detail, we wish to replace the valence bonds jl^jj^x CT2 by bosons Bj air72 , which 
fulfill exact bosonlike commutation relations: 



R R^ 

- D i,CTl<T2 > JD j,<J3<J4 



^ (71, (73 < ^0"2,CT4 • (4) 



There are many ways to implement such a replacement. In the Generalized Holstein- 
Primakoff (GHP) approach ||, which we follow here, the mapping is defined by imposing 
the requirement that all quadratic operators in the original space are mapped in such a way 
as to preserve their commutation relations. More specifically: 

• The boson image of quadratic operators are assumed to be given by Taylor- series 
expansions: 

F B = + + + (5) 



The terms in these Taylor expansions are obtained via the condition that any com- 
mutation rule [A, B] = C between the original set of quadratic operators must be 
preserved at each order of the expansion: 



A<®,B<»] + \aw,b<® 



c (o) 



(6) 



We should note here that boson mappings have typically been applied to systems of 
interacting fermions. However, there is no fundamental difficulty in applying them to systems 
of interacting bosons [fLOfl , as in the problems under discussion. 

We have succeeded in building an appropriate boson mapping of the Holstein-Primakoff 
type for valence bonds. In the present discussion, we simply present the relevant mapping 
equations, leaving more detailed discussion of the formalism to a subsequent publication. 

The full algebra of quadratic operators in the Schwinger boson space includes both 
particle-hole (p-h) operators of the form % a ^j a ' and particle-particle operators of the form 
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For the purpose of treating the dynamics of spin-chain hamiltonians, we only need to 
know how to map the on-site p-h operator '-fj and the p-p operators involving bonds 
between neighboring sites, ij^ll+iy an d li,ali+i,a' ■ The relevant images through first-order 
in the Taylor series expansion are as follows: 

The boson image of the on-site particle-hole operator is: 

{l\,a\li,^2) B — ^ {-^,0103-^,0203 + B\_ l (T3ITl 0-30-2 } • (7) 

0-3 

Note that the on-site p-h operator maps exactly onto a p-h operator in the ideal boson space, 
with no need for a series expansion. 

The p-p operators involving nearest neighbor bonds do require infinite series expansions. 
The lowest (zeroth) order images are straightforwardly given by 

f 7 f 7 t y o) =J B? (8) 

y H,<Tl (1+1,(72 J £j l,0l02 \ J 

and 

The first-order images are more complex, as they provide the first reflection of exchange 
effects of the (Schwinger boson) constituents between neighbor bonds. The required first- 
order image that fulfills the second line of (6) is given by 

(^7i,oi7i+l,02 J B ~ 7} 2.^1 {^,0-103-^1,0401 -^,0403 + -^,0103-^+1,0204-^+1,0304 
03,04 

+ jB l,0302 5 l-l,0 4 01 ^-1,0403} • (10) 

The corresponding image associated with annihilation of a bond between nearest neighbors 
is obtained from (10) by hermitian conjugation. 

In fact, closing the algebra to first order requires inclusion of the p-h operator between 
next-to-nearest neighbor [i with i + 2) sites. We will not discuss this any further here, as it 
does not impact on the analysis or results to follow. 

There is one further complication that should be noted before considering the appli- 
cation of these methods. The mapping equations given above can be applied in several 
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different ways to a given spin-chain hamiltonian. One possibility is to express the two-body 
interaction entering the hamiltonian in p-p form (7^7^77) and to map using eqs. (8-10). 
Alternatively, the hamiltonian could first be transformed into p-h form (7^77^7) and then 
mapped with eq. (7). A third possibility, of course, is to map part of the hamiltonian in 
p-h form and part in p-p form. Were we to do the resulting analysis exactly in the ideal 
boson space, all such approaches would be equivalent. In variational treatments, on the 
other hand, it is essential to map the hamiltonian in such a way as to maintain the key 
correlation effects. 

As a specific application of these methods, we consider a spin-chain system governed by 
the alternating bond hamiltonian 



This system has been studied extensively, especially with regards to its dimer phase transi- 
tions. 

The variational treatment we will apply to this system is based on singlet bonds only. 
Thus, following our earlier remarks, we must map the hamiltonian so as to most efficiently 
reflect these bonds. The appropriate separation is: 



N 




(11) 



i=l 




(12) 



where 



N 



Hi(a) = £[1 - {-ya]S?St 



i=i 




(13) 



and 





(14) 
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We map the term Hi(a), corresponding to SfS? +1 , in p-h form and the term H 2 (a), corre- 
sponding to S^Sf +1 + SfS^, in p-p form. 

Applying the mapping in this way to the alternating bond hamiltonian (11) and then 
projecting onto singlet bosons, defined by 



a] 



^[B{ + _-Bl_ + \ , (15) 



we obtain 
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1 N 

H B (a) = - jYll 1 ~ (-T a ] faUi + + + crl-i^-iki} • (16) 

4 i=i 

The trial wave function we use in our variational description of the alternating bond 
spin-chain system is 

N t n a t "e 

l$n ,n e > = II \ ,^ |0>, (17) 
i(odd)=l V n o- n e- 

subject to the Schwinger constraint 

n + n e = 2S . (18) 

This trial wave function reflects the various phases of the system. The Heisenberg phase 
corresponds to n D = n e , with all sites involved in an equal number of bonds with its near- 
est neighbors on each side. The corresponding state is translationally invariant and is the 
analogue of the VBS state in the ideal boson space. Increasing n Q or equivalently n e corre- 
sponds to successive partial dimerization. The cases in which either n Q = or n e = involve 
complete dimerization. 

Here we focus our analysis on the location of the critical points associated with a tran- 
sition from one phase to another. Defining 

£n ,n e (a) = < ®n ,n e \H B (a)\$ nojne > , (19) 

the critical points are given by the condition 

£n ,n e = £n ±l,7i e =Fl ■ (20) 



Straightforward analysis yields for the energy functional 

N 



and for the critical values of a, 



1 + a){— H Yn n e \ 

+ (1 - a )i— + g n ° n ^ 



(21) 



« = 2n ; ( ^ 1) ^ - n o = 0,l,...,2S-l. (22) 

In Table 1, we present the results of this analysis for several values of the spin S. We 
compare the results obtained from our simple analytic formula (22) for the crossing points 
with those from "exact" calculations fTT|-|i~3|| and from the NLSM [|1],|3| . Our simple formula, 
obtained by using a first-order boson mapping treatment and a simple product trial state, 
yields perfect agreement with the exact results where available. This is to be contrasted with 
the NLSM results, which are in much worse agreement. Note that we have also included in 
the table predictions for the location of the phase transition points at higher spins. 

The fact that our first-order boson mapping works so well in describing the location of 
dimer phase transitions does not mean that higher-order contributions are of no dynamical 
importance. Indeed, at the same level of approximation, the energy per site of the pure S = 1 
Heisenberg lattice (a = 0) is calculated to be —5/4, whereas the exact result is —1.401... ||. 
Thus, there is clearly room for improvement, which we expect would be provided in part by 
the next-order contribution of the boson mapping and in part by fluctuation effects. What 
is important to realize, however, is that the method we have outlined provides a systematic 
way to incorporate such improvements. 

There are several areas in which we expect these methods to be useful in the future. Our 
immediate plan is to generalize the hamiltonian (11) to include crystal fields and a uniform 
magnetic field and to study the phase diagram, excitations, magnetization curves, etc. 

This work was supported in part by the DIGICYT (Spain) under contract no. PB95/0123 
and by the National Science Foundation under grant # PHY-9600445. Fruitful discussions 
with Siu Tat Chui and German Sierra are gratefully acknowledged. 
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TABLES 

TABLE I. Location of the crossing points for the alternating bond spin chain. In addition to 
the results of the present analysis, we present results from the non-linear sigma model (NLSM) and 
from "exact" numerical solution. In the case of the exact analyses, only the nonnegative crossing 
points are shown. 



S Exact Present Results NLSM 

1 0.25 ±0.01 ±1/4 ±1/2 
3/2 0. , 0.42 ±0.02 , ±2/5 , ±2/3 

2 0.05< a <0.3 , 0.5< a <0.6 ±1/6 , ±1/2 ±1/4 , ±3/4 
5/2 , ±1/3 , ±2/3 , ±2/5 , ±4/5 

3 ±1/7 , ±3/7 , ±5/7 ±1/6 , ±1/2 , ±5/6 



a The exact results for S = 1 are from Ref. 11, those for ,5 = 3/2 are from Ref. 12, and those for 
S = 2 are from Ref. 13 
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